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Abstract. A framework is a graph and a map from its vertices to R d . A framework is called 
universally rigid if there is no other framework with the same graph and edge lengths in R d for any 
d' . A framework attachment is a framework constructed by joining two frameworks on a subset of 
vertices. We consider an attachment of two universally rigid frameworks that are in general position 
in R d . We show that the number of vertices in the overlap between the two frameworks must be 
sufficiently large in order for the attachment to remain universally rigid. Furthermore, it is shown 
that universal rigidity of such frameworks is preserved even after removing certain edges. Given 
positive semidefinite stress matrices for each of the two initial frameworks, we analytically derive 
the PSD stress matrices for the combined and edge-reduced frameworks. One of the benefits of the 
results is that they provide a general method for generating new universally rigid frameworks. 

1. Introduction 

A framework in M. d is a graph embedded in Euclidean (i-dimensional space. The graph vertices 
are assigned point coordinates, the edges are represented as line segments connecting corresponding 
points, and the edge lengths are determined as the distances between the edge-connected points 
in M d . One of the important questions in geometry of frameworks is whether there are other 
frameworks with the same edge lengths. When searching for such frameworks we want to exclude 
shape-preserving trivial euclidean transformations such as translations, rotations and reflections, 
applied to already found frameworks. Thus different frameworks will differ in pairwise distances 
between some pairs of vertices that are not connected by an edge. In case there are no different 
frameworks in Mr, the original framework is called globally rigid. Furthermore, if there are no 
different frameworks in any space of higher dimension, the framework is called universally rigid. 

Connelly first studied universal rigidity in the context of Cauchy polygons [3]. Later, in [5], 
he set and proved a sufficient condition for an arbitrary framework to be universally rigid. The 
condition is formulated in terms of stresses acting on the edges, and it states that a framework in 
M. d is universally rigid if there exists a positive semidefinite stress matrix of nullity d + 1 (see below 
for definitions). Recently, Gortler and Thurston [8] proved that this condition is also necessary for 
generic frameworks, in which the vertices' coordinates are algebraically independent over rationals. 
Alfakih [9] then established sufficiency of the condition for the wider class of frameworks that are 
in general position, such that no subset of d + 1 vertices is affinely dependent. In [8], Gortler 
and Thurston have also suggested an algorithm for searching for such PSD stress matrices based 
on semidefinite programming. (See [TB], [17] for related work, and [13] for a general overview of 
semidefinite programming.) 

Although universally rigid frameworks are now characterized by the above conditions, not many 
results are known about operations on frameworks that preserve universal rigidity, nor are there 
many known examples of systematically generated universally rigid frameworks. In [3], Connelly 
generated larger Cauchy polygons inductively from smaller ones. Another known framework exten- 
sion method that preserves universal rigidity is (d+ l)-lateration |15| jlOj, which is also too specific 
and not efficient in terms of the number of edges added to the original framework. 



I would like to thank Dylan Thurston for his tremendous help and support, and Andrew Fanoe for fruitful 
discussions. This work was partially supported by NSF Grant DMS-0739392. 

1 



2 



KIRIL RATMANSKI 




Figure 1. Attachment of two universally rigid frameworks in M 3 . Vertices 1,2,3 
and 4 are shared by the two frameworks. Compared to Ga, graph Gb has an 



additional edge {1,3} between the shared vertices. By Theorem 15, removing this 
edge from the attachment G preserves its universal rigidity. 

In this paper we consider a new problem of characterizing the universal rigidity of attachments 
of two universally rigid frameworks in general position in W 1 . By attachment we mean a combined 
framework in which the two frameworks share some vertices and in which all edges are preserved 
(see Figure 0. We prove that an attachment is universally rigid if and only if the number of 



shared vertices (at which the two frameworks attach) is d + 1 or greater (see Theorem 12). We also 
show that removing those edges connecting shared vertices that are inherited from only one of the 
attached frameworks preserves universal rigidity of the attachment. Finally, we derive PSD stress 
matrices of nullity d + 1 for both original and edge-reduced attachments explicitly, from the PSD 
stress matrices that are given for the two joined frameworks. 

Thus, similar to frameworks, we can characterize the framework attachments by the stress matri- 
ces. Moreover, the stress matrices for universally rigid framework attachments have the same prop- 
erties as the stress matrices for universally rigid frameworks, and computing such matrices does not 
involve computational search. These results provide a more general method than (d + l)-lateration 
|10| for generating new universally rigid frameworks from two or more arbitrary universally rigid 
frameworks. It can also be used to analyze complex frameworks by decomposing them into smaller 
attached frameworks for which universal rigidity is known (or can be more easily established). The 
results may have a potential application to structural engineering and molecular biology. 

The rest of the paper is organized as following. Definitions and preliminaries are given in section 
[2] In section |3] we formulate and prove the main result which is a necessary and sufficient condition 
for an attachment to be universally rigid. We also prove universal rigidity of the edge-reduced 
framework attachment. Section [4] contains derivation of PSD stress matrices for each of the two 
attachments. In Section [5] we show that edge-reduced framework attachment is a generalization of 
(d + l)-lateration as a method for creating new universally rigid frameworks possessing PSD stress 
matrices of nullity d+1. We also show that an edge-reduced graph attachment of two graphs whose 
frameworks are universally rigid in any general position always has a a PSD matrix of nullity d+1. 



2. Background 

A graph G = (V, E) is a finite set of vertices V = (1,2, ... ,v) together with a set of edges E. 
The set E can be represented as a set of two-element subsets {i,j} of V. We denote the number 
of graph vertices and edges by v = \ V\ and e = \E\, respectively. A graph is called k- connected if 
removing at least k vertices disconnects the remaining vertices. Having two graphs G a = (Va, Ea) 
and Gb = (Vb,Eb) such that Va H Vb 7^ after appropriate labeling of vertices by integers, we 
can construct a graph attachment G = (V,E) where V = Va*JVb and E = Ea U Eb (Figure [T]). 
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Throughout the paper, we assume that the set Va H Vb of shared vertices is a proper subset of both 
Va and Vb- 

A configuration of graph G = (V,E) in M. d is a mapping of V to M. d . A configuration can also be 
represented as a single point p = (pi, . . . , p„) G M rf x • • • x = W" d , where pi is the coordinate 
of vertex i in M. d . A configuration is called generic if its coordinates (in W" d ) are algebraically 
independent over Q. A configuration is said to be in general position if every subset {p^, . . . , Pi d+1 } 
of d + 1 vertices is afnnely independent. Two configurations p and q are congruent, or p = q, if 
q = Rp, where R is an element of the group Euc(d) of rigid motions in M. d including translations 
and rotations. A graph G together with its configuration p is called framework in M d , denoted as 
G(p). Two frameworks G(p) and G(q) are said to be equivalent if whenever is an edge of G, 

then ||pj — pj|| = ||qj — qj||. A framework attachment is a framework of graph attachment defined 
above. Alternatively, if we are given two frameworks both having subsets of vertices with the same 
pairwise distances, it is possible to coincide the two subsets of vertices in space by applying a rigid 
motion to one of the frameworks. A framework attachment is created by merging some or all pairs 
of the coinciding vertices, such that all edges from both frameworks are preserved. There may 
be more than one attachment constructed from the two frameworks. While the remainder of this 
section covers single graphs and related notions, we will treat framework attachments in detail in 
the next sections. 

A framework G(p) in R d is called globally rigid if whenever G(q) is equivalent to G(p) then 
p = q. A weaker than global rigidity is the notion of local rigidity. A framework G(p) in R d is 
called locally rigid if there is e > such that whenever G(q) is equivalent to G(p) and q E B e (p), 
then p = q. In other words, locally rigid framework is globally rigid in some e-neighborhood. By 
definition, global rigidity implies local rigidity. 

A framework G(p) in M d is called universally rigid if whenever G(q) in M rf is equivalent to G(p) 
and d' > d, then p = q. A graph is called universally rigid in ~R d if any of its frameworks in general 
position M. d is universally rigid. Universal rigidity implies global rigidity. 

A necessary condition for global rigidity of a framework is given by the following variation of 
the Hendrickson's theorem [2], which holds for frameworks in general position rather than generic 
frameworks. 

Theorem 1. A globally rigid framework in general position in WL d is (d + V)-connected. 

Proof. See [2], where the same argument for generic frameworks applies to frameworks in general 
position as well. □ 

For a given graph G = (V, E) with v vertices, e edges and dimension d, we define the edge 
function f : W d — > M e as a function mapping graph configuration in W d to its edge-square-length 
in G(p) (after some ordering of e edges): f(p) = /(pi, . . . , p v ) = (..., |||pt — Pj || 2 , . . •)• The exvd 
Jacobian matrix corresponding to the edge function is called rigidity matrix, denoted by df: 

df=^0 ■■■ Pi-Pj ••• Pj -p 4 ••• oj. (2.1) 

Given the rigidity matrix df, a framework is infinitesimally rigid if rank(c|f) = vd — C^ 1 )- Equiv- 
alently, infinitesimal rigidity of a framework can be characterized by its infinitesimal flexes. A 
vector q 6 M. vd is called an infinitesimal flex of the framework G(p) if for any edge {i,j}, the scalar 
product (pi — pj , qj — cy) = 0. A vector q is called a trivial infinitesimal flex of the framework G(p) 
if there exists a differentiable path h(t) of rigid motions in M d such that h(0) = p and h'(0) = q. 
Along the path all pairwise distances (and edge lengths in particular) remain constant, therefore 
evaluating the derivative of (hj(i) — hj(t)) 2 at t = shows that such q satisfies the scalar product 



-1 



KIRIL RATMANSKI 



equation above. A framework G(p) is then called infinitesimally rigid if every infinitesimal flex of 
G(p) is trivial. 

The next notion is key for the main results of the paper. For each edge {i,j} of graph G we 
define a scalar wij = Wji, and have them arranged into a vector w = (. . . , w^, . . .) G M e . We 
say that vector w is an equilibrium stress (or simply stress in this paper), if the following vector 
equation holds for each vertex i: 

J2 ?%(Pi-Pi) = 0. (2.2) 

By writing the above as a system of vd scalar equations and then forming into a matrix form, we 
arrive at the following relation between stresses and the rigidity matrix: 

Proposition 2. The space of stresses equals to k.er(df T ). 

For each element w = (. . . , w^, . . .) of the space of stresses, there corresponds a v x v symmetric 
stress matrix defined such that for i,j G V and i ^ j, Vtij = Wij for {i,j} E E and = 
for {i,j} ^ E. The diagonal entries are defined such that the rows and columns sum to zero: 
= An equivalent definition of the stress matrix, which will be used here, is the 

following. 

Definition 3. A stress matrix of a framework is a v x v matrix Q such that: 

(1) for i,j G V, Ojj = Cljj; 

(2) for i,j G V, i / j and {i,j} $ E, Qy = 0; 

(3) E ie y Oij = for all i G V; 

( 4 ) E je y • Pi = for a11 * G y ; 

Properties (l)-(3) follow directly from the construction of stress matrix described above. In 

addition, ]T\ Pj = Y^j^i 'Pj + ' Pi = Ysj^i Pj- Y^j^i 'Pi = 12 w v (Pj _ Pi) = °> 
which verifies (4). 

The following two theorems establish characterization of universal rigidity by the stress matrix. 

Theorem 4 (Connelly [5], Gortler-Thurston [S]). A generic framework G(p) inM. d , having at least 
d + 2 vertices, is universally rigid if and only if it has a positive semidefinite stress matrix with 
nullity d+1. 

Theorem 5 (Alfakih |9j). Let G be a framework in general position in M. d , having at least d + 2 
vertices. If there exists a positive semidefinite stress matrix of G of nullity d+1, then G is 
universally rigid. 

We will conclude this section by developing some properties of the stress matrices that will be 
needed later. 

Lemma 6. Let A be a d x v matrix having pi,. . . ,p v as its columns. Assume a matrix B is 
obtained from A by appending a row (1, . . . , 1). Then p is in general position if and only if every 
(d+1) x (d + 1) submatrix of B has full rank. 

Proof. If such a submatrix, having columns Pi k , k = 1, . . . , d+1, is not full rank, then ^ ak 

for some non-trivial {ofe}- But this holds if and only if 12 a kPi k = ari d 12 a k = 0) i- e - when 
Pi k are affinely dependent. □ 

Corollary 7. The set of configurations in general position is an open set. 

Lemma 8. If Q is a PSD stress matrix of nullity corresponding to a framework G(p) in 

general position, then the coordinate projections of p and the vector (1, . . . , 1) T are a basis for 
ker(O). 
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Proof. We note that the equation in (4) of Definition^ is a vector equation, and so it is satisfied for 
each of the coordinate projections. Vector (1, . . . , 1) satisfies the same equation by (3). Taking 
d coordinate projections' column vectors, vector (1, . . . , 1) T and forming into a matrix, we get a 
matrix B (as in Lemma [6| transposed, which is of full rank. Therefore the coordinate projections 
and vector (1, . . . , 1) T are d + 1 linearly independent vectors in ker(O), and hence they are the 
basis. □ 

In the following lemma we assume a framework in general position has a PSD stress matrix of 
nullity d+1. 

Lemma 9. For a framework in general position, its PSD stress matrix of nullity d+1 can be 
diagonalized by an invertible matrix whose first d+1 columns are coordinate projections of p and 
a vector (1, . . . , 1) T . 

Proof. Let f2 be the matrix. Since Q is symmetric, f2 = UAU T or, Q.U = UA, where U is a 
unitary matrix of eigenvectors corresponding to eigenvalues of Q arranged on the diagonal of A 
Because O is a PSD matrix of nullity d + 1, A has d + 1 zeros and v — d — 1 positive values on 
its diagonal. By re-ordering columns of U and corresponding diagonal entries of A we make the first 
d+1 columns of U to be the eigenvectors corresponding to O-eigenvalues which are the first d+1 
entries of the diagonal of A. We then replace these d+1 columns of U with (non-orthonormal) 
set of d coordinate projection vectors and a vector (1, . . . , 1) T , which are also the eigenvectors 
corresponding to O-eigenvalues. This results in a new matrix S such that £IS = SA. Since the first 
d+1 columns of U and S are a basis for the same subspace (in particular, each of the first d + 1 
columns of S can be expressed as a linear combination of the first d+1 columns of U) , it is easy 
to show that the set of all columns of S is linearly independent, which implies that S is invertible. 
Therefore O = SAS' 1 . □ 

For general v x v positive semidefinite matrices the following is also true: 

Lemma 10. Let A, B be PSD matrices. Then ker(A + B) = ker(A) n ker(B). 

Proof. If x € ker(A)nker(£), then (A+B)x = Ax+Bx = 0+0 = 0, so ker(A)nker(5) C ker(A+B). 
Conversely, if x € ker(^4 + B), = x T (A + B)x = x T Ax + x T Bx, therefore x T Ax = x T Bx = 0. 
Since x T Ax = {A l / 2 x, A x l 2 x), where A 1 / 2 = UA l / 2 U T and hence ker(A 1 / 2 ) = ker(A), we have 
= A l / 2 x = Ax, and so x G ker(^4). Similarly, x G ker(B), hence kei{A+ B) C ker(j4)nker(i3). □ 

Lemma 11. For any PSD matrix Q\ and a matrix f?2 satisfying ker(Qi) C ker(02), there is c > 
such that c£li + is a PSD matrix with the same nullity as 

Proof. For a sufficiently small e, we can view f2i + e^2 as perturbed matrix Q\. Assume that f^i 
has rank r. Then, the matrix has r positive eigenvalues A; [TT], while the rest of the eigenvalues are 
0. We can assume the eigenvalues are ordered such that A/ > A/ + i > for I = 1 . . . 1 — 1. By Weyl 
theorem [T3], the perturbation of the eigenvalues of f^i is bounded: |A; — A^| < e||^2|| for all I. A 
sufficient condition for preserving positivity of the first r perturbed eigenvalues is |Aj — AJ| < Ar r i, 

where A[ r ] is minimal eigenvalue among the first r. Therefore by setting e = ^[ifM ' ^ e ^ rs ^ r 
eigenvalues of + ef^ are positive. 

The rest of the eigenvalues of f2i + 6^2 are zero. If this was not the case, then rank(S7i + ef^) > 
rank(J)i), or, dim(ker(f2i + ef^)) < dim(ker(Oi)). From other side, since ker(f2i) C ker(r?2) = 
ker(e^2) for e > 0, we must have ker(Oi) C ker(f2i + eOg), which implies dim(ker(Oi + ef^)) > 
dim(ker(f2i)), a contradiction. Therefore Cli + e^2 has r positive eigenvalues while the rest are 
zero, and so the matrix is positive semidefinite with the same nullity as Q\. By scaling the matrix 
by c = e _1 , the conclusion of the lemma follows. □ 
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3. Framework attachments: main results 

In this section we will address the main question of this paper: given two universally rigid 
frameworks, under which conditions their attachment is also universally rigid? 

Theorem 12. A framework attachment of two universally rigid frameworks in general position in 
M. d , not joined on all vertices, is universally rigid if and only if the number of shared vertices is 
greater or equal d + 1 . 

Proof. Assume that we are given the attachment G in M. d of two universally rigid frameworks, Ga 
and Gb, that are in general position. Denote the set of vertices in Ga by Va, the set of vertices 
of Gb by Vb, and the set of shared vertices by Va = Va ri Vb- If \Vs\ < d + 1, removing the 
shared vertices from G leaves a disconnected graph, which means G is not (d + reconnected and, 
by Theorem [TJ not universally rigid. 

To prove the converse, we express can express the configuration p of G as p = (pA)Pb\c)> 
where p^ is the configuration of Ga- Since Ga is universally rigid, any configuration p' of G 
satisfies p^ = p^. Therefore for any such p' there corresponds (by a rigid motion) a configuration 
p" = (j>Ai p'b\c) suc b that p" = p'. We will now show that any such p" is congruent to p, thus 
proving p' = p. 

Without loss of generality, assume that the two frameworks are joined at n > d + 1 vertices 
vi,. . . ,v n . Since each framework is in general position, the coordinate vectors pi, . . . , Pd+i are 
afnnely independent, i.e. yi = P2 — Pi, • • • , yd = Pd+i — Pi are linearly independent. Since Gb 
is universally rigid, p B = p_b- That means there is a congruency map H : ps — > p" B . The map 
fixes pi, . . . , Pd+i, therefore it must be some linear transformation A from the orthogonal group Od 
(we can assume, without loss of generality, that vertex v\ is at the origin). The transformation A 
also fixes the basis vectors yi, . . . , y^. Therefore, for a coordinate vector pj = ^j^iYi °f a vertex 
Vj £ Vb\Ci we have: 

Pj = Apj = a iYi = a i( A Vi) = ^2 = Pi' 

i i i 

which shows that p" B \n = Pb\c- Therefore, p" = p and hence p' = p. Since p' is an arbitrary 
configuration of the attachment, and the above argument did not depend on the dimension d, this 
shows that the attachment is universally rigid. □ 

With this result, one may naturally ask: 

Question 13. What are the generating graphs for universally rigid graphs under graph attachment? 

In other words, we are interested in identifying a smallest set of universally rigid graphs which 
would span, under finite attachment, the space of universally rigid graph attachments. At a first 
glance the space may resemble a semigroup with a set of generators. However, since there can be 
more than one way to attach a pair of graphs, the binary operation is not well defined. Leaving the 
attempt of finding the answer to the above question for a separate possible research, we will now 
establish the second main result, which states that removing certain edges between shared vertices 
in the attachment preserves universal rigidity (see last diagram in Figure [TJ. 

Definition 14. An edge-reduced framework attachment is a framework obtained from framework 
attachment by removing those edges between the shared vertices that are inherited from only one 
of the two joined frameworks. 

Theorem 15. An edge-reduced framework attachment of two universally rigid frameworks is uni- 
versally rigid. 

Proof. Given the non edge-reduced attachment G in M. d and the corresponding edge function /, we 
again denote the set of vertices in Ga by Va, the set of vertices of Gb by Vb, and the set of shared 
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vertices by Vq- As in the proof of Theorem 12 here it will also be sufficient to consider only the 
configurations of G in which the vertices of Ga are "pinned". Let w, defined by p,4 = vr(p), be 
the projection of the attachment configuration p = (pa,Pb\c) onto coordinates of the vertices of 



Ga- Since G is universally rigid (by Theorem 12), the set / 1 (/(p)) of all edge-length preserving 



configurations consists of configurations congruent to p. 

Assume now we remove the edge {i, j} of length l\ between shared vertices that was inherited from 
Gb only, resulting in an edge-reduced framework attachment Gi with a corresponding edge function 
fx. Since Ga is universally rigid and the removed edge did not belong to Ga, any configuration 
p' = (p^P^y^) 6 /r 1 (/i(p)) satisfies p' A = pa- Therefore for any such p' there corresponds (by 
a rigid motion) a configuration p" = {pa,p" b \ c ) £ /j" (/i(p)) H it~ 1 (pa) such that p" = p'. We 
will now show that any such p" is congruent to p, thus proving p' = p. 

We make the following observation regarding the configuration sets: 

r 1 (/(p))c/ r i (/ 1 ( P )), 
r 1 (/(p)) = /r 1 (/i(p))nA i , 



where Dij is the set of all edge-reduced configurations with the distance between vertices i and j 
equal to lx- By restricting the attachment configurations to those with vertices of Ga mapped to 
Pa, we get: 

p e r x (/(p)) n vt-^pa) = fxHh(p)) n n n~\p A ) = fr\h(p)) n ^(pa) i p", 

where the last equality stems from the fact that if pa is fixed then dist(i, j) = lx, or tt~ 1 (pa) C D^, 
and so the edge length restriction can be omitted. Therefore p = p" = p', and since p' is an 
arbitrary configuration of Gx, this shows that Gi is universally rigid. Proceeding by induction, 
we conclude that removing those edges between shared vertices that are inherited only from Gb 
preserves universal rigidity. □ 



The main idea here is that universal rigidity of Ga maintains the same pairwise distances between 
the shared vertices for all edge-reduced attachment configurations. Since G is also universally rigid, 
we conclude that all edge-reduced attachments agree on pairwise distances between all vertices, and 
so, in particular, all configurations are congruent to p. 



4. Stress matrices for framework attachments 

Let's assume we are given an attachment of two universally rigid frameworks in general position 
with PSD stress matrices Qa of size va xva, and £Ib of size vb x vb, each of of nullity d + \. Here 
va and vb are the number of vertices for the corresponding frameworks, and so the framework 
attachment will have v = va + vb — n vertices. 

We re-order, if needed, the columns and rows of £Ia such that the last n columns (and rows) 
correspond to shared vertices. We also rearrange fi^ to have the first n columns (and rows) 
correspond to shared vertices. Next, we extend framework Ga by vs — n non-shared vertices of Gb 
(leaving out new edges), thus getting an extended framework Ga- Adding the new disconnected 
vertices is equivalent to augmenting £Ia by vb — n rows and columns of zeros stresses, which results 
in an extended v xv stress matrix CI a- Similarly, we extend framework Gb to Gb by adding va — n 
non-shared vertices (and no edges) of framework Ga- The stress matrix CIb is extended to a v x v 
matrix fig in a similar fashion except for the way we augment the matrix by zeros, as shown below. 
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Figure 2. A framework in R 2 of a Tutte graph |4|. The framework is universally 
rigid and has a PSD stress matrix of nullity 3, but it is not infinitesimally rigid. The 
arrows indicate a non-trivial infinitesimal flex. The numbers by the edges are the 
stresses corresponding to a PSD stress matrix of nullity 3. 



(The difference comes from maintaining consistent labeling of vertices.) 
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Clearly, the attachment of the two extended frameworks (joined at all vertices) is identical to the 
original attachment. Also, the two extended stress matrices are of the same size. We can now state 
the two main results of this section. 

Theorem 16. Given PSD stress matrices £l A and £Ib of nullity d + 1 for two frameworks in 
general position sharing n > d + 1 vertices, a PSD stress matrix of nullity d + 1 for the framework 
attachment G can be obtained by summing the two matrices after extending each by appropriate 
number of zero columns and rows: 



( 



\ 



o 



o 



\ 



+ 



n 



B 



(4.1) 



/ 



Theorem 17. IfG A is also infinitesimally rigid, a positive semidefinite stress matrix of nullity d+1 
for the edge-reduced framework attachment, in which K edges inherited from Gb were removed, is 
of the form: 

Q re = cVt A + VL A k + Q B , 

where c is a large enough constant, and the entries of matrix £1 A k o,re derived from the rigidity 
matrix of G A and stress matrix by solving K systems of linear equations. 

In the rest of this section we prove these two results. 

Remark 18. While being complementary to the two main theorems from the previous section, these 
results do not actually require universal rigidity of the attached frameworks. From other side, a 
condition of infinitesimal rigidity of the framework, whose edges are preserved in the edge-reduced 
attachment, has been introduced for the construction of the stress matrix in the proof of Theorem 
17 An example in Figure [2] shows that infinitesimal rigidity is not automatically satisfied even 
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when the framework is universally rigid with a PSD stress matrix of nullity d + 1. Infinitesimal 
rigidity does, however, hold for all generic universally rigid frameworks. 



Proof of Theorem 16. The matrix f2 is a PSD matrix because both Qa and £Ib are. Also, since 
both Qa and Qb satisfy properties (1), (3) and (4) of Definition |3j so does ft. If i,j 6 Va U Vb, 
i / j and {i, j} ^ Ea U -Kb then the entry of both Qa and Qb is zero, and so is entry 

of f2. This verifies property (2) as well, and hence the matrix is a PSD stress matrix for the 
framework attachment. 

To prove that f2 is of nullity d+1, we first consider the case where the number of shared vertices 
is exactly d+1. Assume, as above, that the last d + 1 columns(rows) of £Ia, and the first d + 1 
columns (rows) of £Ib correspond to the shared vertices, in the same order. We begin by diagonaliz- 
ing the PSD matrices by Ua = UaAaU^ and Qb = SbAbS^ 1 , where £Ia and £Ib are diagonal, Ua 
is unitary, and Sb is the invertible matrix constructed as in Lemma[9| The corresponding extended 
stress matrices, obtained as in the same lemma, are then readily diagonalizable as £Ia = UaAaUa 
and &b = Sb^-bS^ 1 , where the v x v matrices are of the form: 



A A 



U a 
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Sb 














J 



Using the above diagonalization, we can write: 

n = n A + h B = u A (A A + u^bUaM = u a (a a + v^SbAbS^UaWL 



or, 



n = U A (A A + v- l l B v)ul, 

where V = S^Ua- Since Ua is invertible, the matrices (A^ + V~ 1 AbV) and O are of the same 

nullity (fi) = dim(ker(A J 4 + 
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nullity. Also, both A^ and V 1 AbV are PSD, and so, by Lemma 
V~ 1 AbV)) = dim(ker(A,4) n ker(U _1 A#y)). Using the formula for the dimension of the sum of 
two vector spaces [12] . we expand the last expression for the intersection, which gives: 

nullity(^) = dim(ker(A4)) + dim(ker(U~ 1 A B U)) - dim(ker(A A ) + ker(U^ 1 A B U)). (4.2) 

We will now determine each of the three terms appearing in the above equation. Since only the 
last vb diagonal entries of A^ are zero, the basis for ker(Ayi) consist of v x 1 vectors ej, where 
i = {va — d,...,v} and ei = (0, . . . , 0, 1, 0, . . . , 0) T , with 1 at the i th coordinate. In a compact form, 



I o N 



ker(A^) = im 



(4.3) 



\ / 



vxvb 
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Similarly, 



(t \ 



V A 



ker(A^) = im 

V J 

VXVa 

and since dim(ker(V _1 A£V)) = dim(ker(As)), we get: 



(4.4) 



dim(ker(A j4 )) = v B , 
dim(ker(F- 1 A B F)) = v A . 



(4.5) 
(4.6) 



For the third term in (4.2), we will find the basis of ker(A^) + ker(V _1 A£V). We first note that 
x £ ker(A^) V~ 1 x 6 ker(V _1 A#y), and since V is non-singular, the basis B of ker(y~ 1 A^y ) 

consist of the columns of the matrix-product of V~ l and the basis-matrix of ker(A#), given in (4.4): 



(i \ 

J-V A 



B = V 



(t \ 



uls B 



1 <-'A 



V / 



V 











Sb 



\ 



(i \ 



\ 1 



By splitting Sb and then U\ into blocks of appropriate sizes to perform block matrix multiplication, 
we get: 



B 



( 


I 


! o I 

-i -i 





\ 




\ 


/ 


ul 










: D \ 

-\ ^ 


Y 



















\x\ 


Z 














I 


V 








) 


V 


/ 


V 







\ 



( 



J 



I 







\ 



\D 



\ 



X 



J 




(4.7) 



where the matrix D is the (d + 1) x (d + 1) upper left block of the matrix Sb- By construction, 
the first d + 1 columns of Sb are the d coordinate projections of Gb and a vector (1, . . . , 1) T , 
and so the first d columns of D consist of coordinates (of shared vertices) while the last column 
is the (d + 1) X 1 vector (1, . . . , 1) T . Since Gb is in general position, by Lemma M the matrix D 

is invertible. This in turn implies that the above va x va submatrix C/J ( ^ Dj * s i nver tibl e - 

Therefore, right-mulitplication of B, given in (4.7), by the inverse of the submatrix produces a new 
basis-matrix [12J for ker(y _1 A#V): 



kei(V~ l A B V) = im 




\ 



U A 



( T. \ 



im 



Q 



V / 



(4.8) 



VXVA 
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for some matrix Q. Combining (4.3) and (4.8) now gives: 



' o \ 



ker(A A ) + ker(y _1 A B V0 = im 



L VB 



( L. \ 



+ im 



V / 



•-V A 



im(I v ) 



and hence 



dim(ker(A A ) + ker(y- 1 A B F)) 



va + vb - d - 1. 



(4.9) 



By substituting (4.5), rt4.6| and aA.9« into the dimension formula (4.2), we finally get: 



nullity (fi) = vb + va — (va + vb — d — 1) = d+ 1, 



which completes the proof for the case where the number of shared vertices is d + 1. 

In case the number of shared vertices is n > d + 1, we consider an "intermediate" framework 
attachment which has d + 1 shared vertices and n — d — 1 pairs of vertices that are coinciding in 
space. The attachment with n shared vertices is then constructed by merging each pair of coinciding 
vertices into one (shared) vertex. Such merging in turn corresponds to certain operations on the 
stress matrix of the intermediate attachment, which we will now describe in detail. 



The PSD stress matrix for the intermediate attachment is similar to (4.1): 



n 



Cl C 2 
+ 

+" | 

'- + 



o \ 



+ 



r-2 



n B J 



where a pair of columns c±, C2 (and rows n, T2) correspond to one of the n — d—1 pairs of coinciding 
vertices that are to be merged. The matrix Q is of nullity d + 1 (which, as in the proof above, comes 
from Gb being in general position). We will now obtain a stress matrix for the attachment with 
the two coinciding vertices merged, and show that this matrix is positive semidefinite of nullity 
d+1. 

We first add column C2 to column c%, and row T2 to row r\. Such operation results in a matrix 
W = mR T where R is an elementary matrix, and so tt' is a PSD matrix of nullity d+1. Adding 
the columns and rows corresponds to merging the two vertices. We now exchange the column C2 
(row r%) with the last column (row), resulting in another, positive semidefinite matrix Q" of nullity 
We can view 0" as a matrix fl m bordered with a vector y and a real number a: 



0" 



17. 



\ 



V y 1 

By theorem 4.3.8 in [11], the eigenvalues /ij of £l m are related to the eigenvalues Aj of ft" by the 
following inequality: 

Ai < m < A 2 < lii < ■ ■ ■ < Xv-i < Hv-\ < K- 
Since O" is of nullity d + 1, we have Ai = • • • = A^+i = 0, and hence [i\ = ■ ■ ■ = fj,^ = 0, 
< [J-d+i < Ad+2, and /J.d+2, ■ ■ ■ , l^v-i > 0. This implies that VL m is of nullity d ov d+1. However, 
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Q m is a stress matrix for a framework attachment (with the pair of coinciding vertices merged) , and 
as such its minimal nullity is d+1. This comes from the fact that the configuration corresponding 
to the attachment of two frameworks in general position with at least d+1 vertices each does not 
lie in any affine subspace of M. d [Sj. Repeating the above operations on the n — d — 2 rows/columns 
of Q m corresponding to the rest of coinciding vertices, we then obtain a PSD stress matrix of nullity 
d+1 for the attachment with n shared vertices. This completes the proof of Theorem fl6| □ 



Before formally proving Theorem 17, we first outline the proof for the case where only one edge 
(say, {i,j}) of Gb is removed from the attachment. We start by having two frameworks and the 
edge {i,j} still connected in Gb- From the PSD stress matrix f2#, we obtain a stress w\ on the 
edge We then add the edge {i, j} to the same pair of vertices in Ga, thus obtaining edge- 

extended framework Ga, and use Proposition |l9| below to derive a stress vector w such that the 



stress on the added edge is —w\ (as in Figure^]). With the help of Lemma 11, we then obtain a 



PSD stress matrix of nullity d+1 bearing the same stress on the added edge. Attaching the two 



frameworks, by Theorem 16, will correspond to summing the corresponding extended PSD stress 



matrices, resulting in w\ + {—w\) = stress on the edge of the framework attachment. By 



Theorem 15, we can remove this edge from the attachment. Since the entry in the stress matrix 



corresponding to {i, j} is 0, after removing this edge the same stress matrix will satisfy equilibrium 



equation 2.2 and hence it is the stress matrix for the edge-reduced attachment as well. In the 



proof of Theorem 17 below, we will follow this outline for a more general case where K edges are 
removed from the attachment. But first we establish the following preliminary result guaranteeing 
a non-zero stress on the added edge. 

Proposition 19. For a framework Ga, obtained from infinitesimally rigid framework Ga by adding 
edge {i,j}, there exists a stress vector w such that the stress on the added edge is equal to —w\. 



Proof. For the unmodified Ga we have the rigidity matrix df, defined by (2.1). By proposition [2j 



a stress vector is an element of a kernel of rigidity matrix transposed. Since Ga is infinitesimally 
rigid, rank(4f) = VAd — f^ 1 )- With being the number of edges (and the number of rows 
in df), nullity {df T ) = — rank(o!f). Adding edge {i,j} to Ga corresponds to adding a row to 
df (or, a column to df T ), resulting in a new rigidity matrix df. Since adding an edge preserves 
infinitesimal rigidity, rank(df) = VAd — ( = rank(d/), and so nuhity(<i/ ) = + 1 — rank(o!f) = 
nullity (o!f T ) + 1. In other words, adding an edge to an infinitesimally rigid framework increases the 
dimension of the space of stresses by 1. We will now show that in the new space of stresses, there 
exists a vector with non-zero stress component corresponding to the added edge We first 

note that columns of df T and components of a stress vector w G ker(df T )) are indexed by existing 
edges of Ga- We append a new edge related column (in df T ) and a component (in w), and from 
the following equalities: 



•V T ™=(jf T \ P ){ w o)=W T i w 

where p = (0, . . . Pi — Pj, . . . , Pj — pj, . . . , 0) T corresponds to the added edge {i, j} (see Figure 3c), 
we have: 

weker(df T ) <=> (o) Gker(df T ). 

Since dim(ker(d/ )) > dim(ker(d/ T )), vectors (w,0) T do not span all of ker(o!f ), and so there 

exists a stress vector w = (w,tp) T E ker(o!f ) such that tp ^ 0. By scaling, we can set ip = —w%, 
while the rest of the components (namely, w) of vector w can be found by solving the equation: 

or, in a compact form, 



dJ T w=( d f p)(_" 1 )=0, 



df 1 w = pwi. (4.10) 
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(A) 



(B) 




\-ULh 



(c) 



Figure 3. Edges between the shared vertices in Ga (3a), Gb (3b) and Ga (3c). 



The stresses on the added edges in Ga counter-balance the stresses on the same 
edges in Gb- 



□ 



Proof of Theorem 11 Following the same construction as in the proof outline above, we first add 



one of the K edges to Ga- From 4.10, we then compute a stress vector for the extended framework, 
whose last entry, corresponding to the stress on the added edge, equals to the sign-inverted stress on 
the same edge in Gg. Repeating for all K edges, we obtain K stress vectors by appropriately 



choosing w\ from Qb and solving 4.10 for each edge. From each Wk, in turn, we form a stress 
matrix such that the (ik,jk) entry, corresponding to the added edge, is equal to the last (non-zero) 
entry of Wk- By summing all K stress matrices, we get a stress matrix VLak for the framework 
Ga constructed by adding all K edges to Ga- The elements of SIak corresponding to the added 
edges will have the same stresses, but with the opposite sign, as the stresses on the same edges 
m. Gb- However, this matrix in general may neither be PSD nor have nullity d + 1. From other 
side, the given matrix O^, which is a stress matrix for both Ga and Ga, is a PSD matrix with 
nullity d + 1, with the entries corresponding to the K added edges equal to 0. Moreover, since 
Qa and SIak both satisfy properties (3) and (4) of Definition [3J while Qa also satisfies Lemma [sj 
we have ker(S7/i) C kei(QAK)- Therefore by Lemma 11, with a suitable choice of c > (such as 



A 



[v A -d-l] 



(Oa)- 



the matrix c£Ia + ^AK is positive semidefinite with nullity d + 1 , and has the 



same stresses on the added edges as Qak- 

Finally, we attach the frameworks Ga and Gb- The K added edges in Ga were already in 
Gb, and so this attachment is identical to the regular framework attachment of Ga and Gb- In 
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terms of matrices, however, this attachment corresponds to adding extended fl^ to the extension 



of cVLa + ^ak ! resulting, as in the proof of Theorem 12, in a PSD stress matrix of nullity d+ 1: 

Its entries, corresponding to the K edges, are equal to 0, therefore Vt re is the desired stress matrix 
for the attachment with the K edges removed. □ 

5. Generation of new universally rigid frameworks 

In ]1Q\, Alfakih et al. proved that (d + l)-lateration framework in general position in R d is not 
only universally rigid but also has a PSD stress matrix of nullity d + 1 . In this section we will show 
that edge- reduced framework attachment is a generalization of (d + l)-lateration, which implies 



that a corresponding PSD stress matrix can be constructed as in Theorem 17 We will also show 
that, in general, the edge-reduced attachment can be used to generate universally rigid graphs. 

Definition 20. A graph G = (V, E) is a (d + l)-lateration graph if: 

(1) verices {1, . . . , d + 1} form a complete graph 

(2) each vertex i > d + 1 connects to d + 1 vertices from the set {1, — 1}. 

The definition above suggests that (d + l)-lateration graphs, after appropriately labeling the 
vertices, can be generated inductively. Starting with a complete graph G, at each induction step 
we expand G by adding a new vertex and connecting it to some d+ 1 existing vertices. The following 
proposition shows that the same can be accomplished using edge-reduced graph attachment. 

Proposition 21. Any (d + l)-lateration graph can be generated by a sequence of edge-reduced 
attachments. 

Proof. Assume at k th step we introduce a new vertex Vk in graph G, and want to connect it to 
d + 1 vertices , . . . , Vi d+1 that are already in G. This can be done by first attaching graph G to a 
complete graph consisting of vertices Vi x , . . . , v j d+1 , Vk such that v i x , . . . , Vi d+1 are shared, and then 
removing the edges between the shared vertices that were not in G. □ 

To obtain a PSD stress matrix of nullity d + 1 for a {d + l)-lateration framework in general 



position we can then use Theorem 17 after checking the condition that any such framework is 
infinitesimally rigid. For this we use the following re-phrased result by Connelly ([1], Proposition 
2.21). 

Proposition 22 (Connelly pfj). Suppose G(p) is an infinitesimally rigid framework in general posi- 
tion in M d . Extend the framework by adding a new vertex Vk and connecting it to d existing vertices 
Vi x , . . . , Vi d such that is not in the affine span of p^ , . . . , Pi d . Then the extended framework is 
infinitesimally rigid. 

From here it immediately follows: 

Corollary 23. Any (d + 1) -lateration framework in general position is infinitesimally rigid. 

With this result we can apply Theorem [17] and compute a PSD stress matrix of nullity d + 1 
for a {d + l)-lateration framework in general position. Since the framework generation involves 
attachment of a complete graph (containing a new vertex) and then removal of some K edges 
between the shared vertices, each step would require computing a PSD stress matrix of nullity d + 1 
for the complete framework, and then solving K systems of linear equations as described in the 
theorem. 

The edge-reduced attachment can also be applied for generation of new universally rigid graphs 
from smaller and simpler graphs whose frameworks in general position are always universally rigid. 
An example of such universally rigid graph, which can not be generated by {d + l)-lateration, is 
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Figure 4. Edge-reduced graph attachment of two universally rigid graphs. Each 
of the two graphs can be generated by 3-lateration, which is not the case for their 
edge-reduced attachment. 

shown in Figure [4j It turns out that all frameworks of the graphs generated by the edge-reduced 
attachment have PSD stress matrices of nullity d + 1, as the following theorem shows. 

Theorem 24. Assume we have two frameworks in general position corresponding to two universally 
rigid graphs. If both frameworks have PSD stress matrices of nullity d + 1, their edge-reduced 
attachment also has one. 



Proof. To apply Theorem 17 we only need to prove that any of the two frameworks is infinitesimally 
rigid. Assume that a framework G(p) in general position of such universally rigid graph G is not 
infinitesimally rigid. Then there exists a non-trivial infinitesimal flex q. By Corollary [7J there 
is t > such that the two frameworks G7(p + tq) and G(p — iq) are in general position. From 
(Pi — Pj > — Qj ) = we have 

UPi + tm) - (pj +tqj)|| 2 = ||pi - Pj|| 2 + \\tqi - tq^ 2 = ||(pj - ten) ~ (Pj ~ tqj)\\ 2 , 

which holds for any edge {i, j}. Therefore, since the flex tq is non-trivial, G(p + tq) and G(p — iq) 
are equivalent but not congruent, which contradicts universal rigidity of the graph G. □ 
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